1 



a^ 

> 
o 

(N 



Gravitational effects on tfie neutrino oscillation in vacuum 

N. Fornengo C. Giunti ^, C. W. Kim ^ and J. Song 

'^Dipartimento di Fisica Teorica, Universita di Torino and INFN, Sezione di Torino, via P. Giuria 1, 
10125 Torino, Italy 

^'Dept. of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland 21218, USA 
and School of Physics, Korea Institute for Advanced Study, Seoul 130-012, Korea 

'^Center for Theoretical Physics, Seoul National University, Seoul 151-742, Korea 

The propagation of neutrinos in a gravitational field is studied by developing a method of calculating a co- 
variant quantum-mechanical phase in a curved space-time. The result is applied to neutrino propagation in the 
Schwarzschild metric. 
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1. Introduction 

The neutrino oscillation phenomenon has been 
discussed in the case of neutrinos propagating in 
a flat space-time both in the plane wave approach 
for highly relativistic neutrinos 0] and in the wave 
packet formalism In this note we present a 
generalization of the plane wave approach to the 
case of curved space-time ■ 

In order to generalize the plane-wave approach 
to the case of a curved space-time, it is useful to 
recall the key points which are used to derive the 
oscillation probability in the flat space-time: 

• Neutrinos are produced at a space-time 
point A{tA,XA) as flavor eigenstates Iva) = 
J2k^akWk)i superpositions of the mass eigen- 
states l^k)- U is the unitary mixing matrix of 
the neutrino fields. 

• The propagation in the space-time of is 
described by the quantum-mechanical phase of 
a plane wave. For a neutrino propagating to 
P{t,x), the state evolves to \i'k{t,x;tA,XA)) = 
exp(— z<f>A;)|i/fc) , where 



= Ek{t- Ia) -pk ■ [x- xa) . 



(1) 



The energy Ek and the momentum pk are related 
by the flat space-time mass-shell relation E'^ — 

Pk^T^l- 



• In order for the oscillation to occur and to be 
observed, the interference among different mass 
eigenstates has to be calculated in a specific 
space-time point B its ,xb)- This gives the space- 
time flavor oscillation probability 

P^^^,.,,{tB,XB;tA,XA) ^^\Uakf\Upk\^+ (2) 

k 

2 ReY,KkUf3kUa,jU;^eM-^^^k3) , 
k>j 



where the phase shift A^kj = ^k 
to the interference between the k^^ and 
eigenstates. 



$j is due 
j**^ mass 



• In actual experiments, the time difference (ts — 
tA) is not measured, whereas the relative posi- 
tion \xB — xa\ of the source and the detector 
is known, and an oscillation probability in space 
P^^-t^^{xB', Xa) can be measured. In the plane 
wave formalism, this can be taken care of consis- 
tently only for relativistic neutrinos by employ- 
ing the light-ray approximation, which consists 
in taking {tB — tA) = \xb —xa\ in the calculation 
of the phase ^k of each mass eigenstate. This cor- 
responds to calculate the quantum-mechanical 
phase of each Ivk) along the light-ray path that 
links A to B. In this approximation, Eq.(|]l) be- 



= {Ek - IPfeDI^B - xa\ 



(3) 



* Presented by N. Fornengo 



Since we are dealing with relativistic neutrinos. 
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(nik Ek), we can approximate, to the first or- 
der, Ek ^ Eo + 0{ml/{2Eo)), where Eq is the 
energy of a massless neutrino. This leads to the 
standard result for the phase shift 



2Eo 



xb ~ xa\ 



2TrLp{A,B) 

T OSC 



(4) 



where Am|^- = to| — rnj, Ll!^^ — (47r£'o/Am^^) is 
the oscillation length and Lp{A,B) = \xb — xa\ 
is the proper length in the flat space-time. 

A generalization to a curved space-time is ob- 
tained by observing that the expression for the 
phase in Eq. (nh is a covariant quantity m : 



where p^'^.^ 



Pi'^dx^ , (5) 



A J A 

— rukdx'^/ds is the four-momentum 
of \vk), g^v is the metric tensor, ds is the line 
element and p^'^' — g^^vP^j,-^ is the conjugate mo- 
mentum to x^ . The last expression in Eq. (^) 
is useful when the metric tensor does not depend 
on some of the coordinates x'^: in that case, the 
corresponding components of p^^^ are constant of 
motion along the classical trajectory of the par- 
ticle, while the components of pj'^j are not con- 
served quantities. The momentum p^^,-^ obeys the 
mass-shell condition in the curved metric 

m^=.g^.p^,,p^,)=.g^>«pW. (6) 

Eq.(^) denotes the quantum-mechanical phase 
acquired by a particle traveling from A to i? in 
a curved space-time. When calculated along the 
classical path, it corresponds to the classical ac- 
tion for a free particle. 

The fact that the time difference between the 
production and detection points is not measured 
leads to the generalization to the curved space- 
time case of the light-ray approximation: the 
phase of each mass eigenstate must be calculated 
along the light -ray trajectory which links Aio B 



(7) 



J light 



This formulation is valid for highly relativistic 
neutrinos, since the light-ray approximation is vi- 
able only in this case. 



2. Neutrino propagation in the Schwarz- 
schild metric 

Let us consider the propagation of neutrinos in 
the Schwarzschild metric. The line element in the 
coordinate frame x^ = {t^r^-d^Lp) is 

ds^ = Bir)At^ — - r^di?^ _ sin^ Mlp^ , (8) 

B{r) 

where B{r) = [1 - (2GA/)/r], G is the Newto- 
nian constant and M is the mass of the source 
of the gravitational field. Since the gravitational 
field is isotropic, the classical orbits are confined 
to a plane, which we choose to be the equatorial 
plane i9 = 7r/2. Therefore, we have di? = 0. The 
relevant components of the canonical momentum 
are p^"^ , p^'' and p^^^ and they are related by the 
mass-shell relation as 

1 / {k)\2 

-^l-^^ipf'^r-BirW^r-^^. (9) 

Since the metric tensor does not depend on the 
coordinates t and tp, the canonical momenta 
Pj'^'* = Ek and plp^ = —Jk are constant of mo- 
tion. With these definitions, the expression of 
the phase is 



dr, (10) 



where Pkir) = —pY"' and the differentials 
(di/dr)o and {d(p/dr)o are calculated along the 
light-ray path from A to B. It is important to 
notice that Ek and Jk, which are constants of 
motion for the geodesic trajectory of l^k), are not 
constant along the light-ray trajectory. Instead, 
the quantities Eq and Jq for a massless particle 
are constant along the light-ray path. 

We briefly comment on the definition of the 
energies involved in the calculation of The 
constant of motion Ek is conveniently used in 
the calculation of the phase but it does not 
directly represent the energy of the neutrino in 
the gravitational field (except for r — oo), since 
Pk{r) = Ek/B{r). Moreover, the correct defini- 
tion of the production and detection energies de- 
pends on the reference frames where the physi- 
cal processes occur. For definiteness, in our dis- 
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cussion we will choose the local reference frame, 
where 4'°'=^(r) = \gur^'^Ek = i3(r)-i/2£;,. 

2.1. Radial propagation 

For neutrinos propagating in the radial direc- 
tion, dip = and (dt/dr)o = ±B{r)^^ where the 
+ (— ) sign refers to outward (inward) propaga- 
tion. Deriving pk{r) from the mass-shell relation 
of Eq.(||) with = 0, the quantum mechanical 
phase of Eq. (nol) becomes 



± 



(^E, - ^ El -B{r)ml)^^. [11) 



We can now apply the relativistic expansion using 
the energy at infinity E^ as a reference value, i.e. 
TOfc < Ek, expanding Eu ~ Eq + 0(to|/(2 £:o)), 
where Eq is the energy at infinity for a massless 
particle. Then, the phase of Eq.(pl]) is easily cal- 
culated and the phase shift which determines the 
oscillation is 



, Ami- 



2Eo 



rA\ 



(12) 



This result does not depend on any assumption 
on the strength of the gravitational field. 

The expression of the phase shift in Eq. ( |l^ ) ap- 
pears identical to that of the flat space-time case, 
Eq.(|^). However, the gravitational effects are im- 
plicitly present in Eq.([T^). When expressed in 
terms of the local energy and the proper distance, 
Eq.(12) becomes (for a weak field, i.e. GM 4C r) 



AmlLp{A,B) 



(13) 



1 - GM 



Lp{A, B) 



where Lp{A, B) ~ — ta + GM ln(r^/rs) is the 
proper distance (in the weak field limit). 

The first parenthesis on the right-hand side in 
Eq.(13) is analogous to the fiat space-time oscil- 
lation phase. The second parenthesis represents 
the correction due to the gravitational effects. 

The proper oscillation length L'^^ which is ob- 
tained from Eq. (|l^) is increased because of the 
gravitational field, as expected. 



2.2. Non radial propagation 

When the classical trajectory is not in the ra- 
dial direction, the motion has an additional an- 
gular dependence. The light-ray differentials are 
(dVdr)o = £;oB(r)-2po(r)^i and (d^/dr)o = 
jQr~'^ B{r)~^PQ{r)^^ . The angular momentum 
is related to E^, the impact parameter h and the 
velocity at infinity 1)^°°^ as Jk = Ekb . Mak- 
ing use of the mass-shell condition Eq.(||) and 
the relativistic expansion ruk Ek, the phase 
becomes 



±- 



'2Eq 



dr 



y/1 - B{r){byr'^) 



(14) 



This expression is valid for any spherically sym- 
metric (and time-independent) field. 

As a specific example, we consider the propa- 
gation of the neutrino around the massive object. 
In the weak field limit, and for the typical situa- 
tion where b <^ ta,Bi the phase shift is 



Ar 



2Ef 



-Lab 



1 - 



2GM 



2rArB Lai 



(15) 



where Lab = {rA + fs)- As discussed for the 
radial case, a comparison of Eq. (^5|) with the flat 
space-time case, requires a correct identification 
of the proper distance and of the energies involved 
in the the physical problem. An interesting ap- 
plication of Eq. to the possibility of gravita- 
tional lensing of neutrinos, has been discussed in 
Ref.i. 
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